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Automated Resummation and Hadron Collider Event Shapes* 

Gavin P. Salam a 

a LPTHE, Univcrsite Pierre et Marie Curie - Paris 6, Universite Denis Diderot - Paris 7, CNRS UMR 
7589, 75252 Paris 75005, France 

This writeup gives an introduction to the theoretical understanding that lies behind automated resummation. 
It then discusses its applications to hadron- collider event shapes. 



> 

o 

<N 
O 
(N 
O 
VO 
O 

22 

Q-h! 

Oh! 
CD 



X 



The all-order resummation of logarithmically 
enhanced terms is essential for the detailed study 
of final-state observables like event shapes and 
jet rates (see for example |T]). The differential 
distribution for an event shape V to have value v 
can be written as 



I da 

a dv 



X'(v)=a s A(v)+a 2 J 2 (v) + 



(1) 



where the leading (LO) and next-to-leading 
(NLO) order coefficients j\ and f 2 can be cal- 
culated with a program such as NLO JET 
At small v, for any observable sensitive to 
soft collinear radiation, the corresponding inte- 
grated distribution E(u) has double-logarithmic 
enhancements at all orders 

2m 1 

E(«)~ Y j Y j H mna ^hx n - ( V ->0) (2a) 

m n— 

= fci(a s L 2 ) + y^r s h 2 (a s L 2 )^ + ... (2b) 

LL S NLLs 

where H mn are numerical coefficients and we have 
introduced the shorthand L = Inl/v. The double 
logarithms are associated with non-cancellation 
between soft and collinear divergences in real and 
virtual graph, since limiting v suppresses real ra- 
diation but not the corresponding virtual terms. 

— 1/2 

When v is small (In v ~ a s ) all terms 
cv™L 2m become of the same magnitude (~ 1) and 
reliable calculations for E(t>) can only be obtained 
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by resumming these terms to all orders, giving the 
leading logarithmic (LL) function h\(a s L 2 ). One 
can systematically reorganise the whole pertur- 
bative series in terms of the logarithmic enhance- 
ments, eq. (I2b|l . with next-to-leading logarithmic 
(NLL) terms, a™L 2m_1 resummed into a func- 
tion yfalh 2 {a s L 2 ), etc., such that each of the h n 
is of order 1. 

— 1/2 

For yet smaller v, L ^> a s , even the reor- 
ganisation of eq. (0) breaks down, since the func- 
tions h n can themselves become large as their ar- 
gument, a s L 2 , grows beyond 1. It turns out how- 
ever that quite often eq. J2J) can be written as the 
exponential of a much less divergent series, 

m+l 

In E(w) ~ G ?nnaTL n (v - 0) (3a) 

m n— 

= L gi (a s L) + g 2 (a s L) +... , (3b) 

LL laS NLL lnT , 

where the key feature of eq. I|3a|l is that the index 
n now runs only up to to + I, i.e. all the double 
logarithms of E arise from the exponentiation of 
a single a s L 2 term in InE. This series can be 
also reorganised systematically, eq. I|3bj) and now 
one redefines the LL terms to be Lgi(a s L) and 
the NLL to be g 2 (a s L), etc., where the g n (a s L) 
have the property that they are all of order I 
as long as L ~ cv^ 1 . This means that the reor- 
ganised perturbative hierarchy remains stable to 
much smaller values of v, than in eq. J2J; further- 
more the NLL terms are more suppressed with 
respect to the LL terms (by 1/L ~ a s ) than was 
the case in the classification of eq. (J2bj). 
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Figure 1. Representation of phase space in e + e — ► 2 jets; emissions (points) in the right (left) hand 
half-plane are collinear to the right (left) hand jet. For further details, see text. 



In this writeup it is the more ambitious, expo- 
nentiated, kind of resummation that will be dis- 
cussed. Let us examine the kinematic plane in 
e + e~, fig. n represented in terms of the trans- 
verse momentum fc t and rapidity (77 = ~ E-k ) °^ 
an emission k, as measured with respect to the 
Born qq axis. The thick yellow lines represent 
the hard-collinear kinematic limit for emissions. 
If there is just a single emission, one calculates 
the O (a s ) term of the logarithmically enhanced 
part of X(i>) by identifying the grey shaded area 
in which the observable's value 2 V{k) is larger 
than v. In the single-emission approximation 
this region is forbidden to real radiation, but 
not to virtual corrections, and the resulting non- 
cancellation leads to a double logarithmic contri- 
bution a s L 2 to E(V). 

In general, when considering multiple emis- 
sions, the boundary V{k\, . . . , k n ) < v will be a 
complex multi-dimensional function of the ki. For 
many observables there are several simplifications 
that can be carried out. These will look quite 
similar to those, based on infrared and collinear 
(IRC) safety, that are used to justify fixed-order 
calculations. We use as a reference scale the lower 
boundary (dashed line) of the shaded region in 
fig. corresponding to a value ~ v = er L of 
the observable for a single emission. Regardless 

2 For compactness we write the observable as a function 
just of the soft and collinear emission momenta, though 
there is also an implicit dependence on the hard momenta. 



of the set of emissions that is present, we require 
that all emissions much softer than the scale of 
the boundary should modify the observable by 
much less than v. These emissions (in grey, be- 
low the dotted line) can therefore be neglected, 
since they will cancel fully against correspond- 
ing virtual terms. Then, only the remaining, 
black emissions, between the dashed and dotted 
lines, will contribute significantly to the observ- 
able. These are clustered in rapidity since there 
is a collinear divergence for the splitting of an 
emitted parton. However, most common observ- 
ables have the property that they are insensitive 
to this collinear splitting, so that each cluster of 
black emissions can be replaced by a single red 
'primary' emission, as shown in fig. ^d. 

For the next stage, it is necessary to exploit 
another property common to most observables, 
namely that 

V(k u . . . , k n ) ~ max[y(*Ei), . . . , V(k n )} , (4) 

so that if V(ki, . . . , k n ) < v then all V(fcj) < v. 
This means that the shaded region will contain no 
emissions (except potentially near its lower edge) 
implying a resummation of purely virtual correc- 
tions there. This leads to the exponentiated dou- 
ble logarithmic contribution Lg\{a s L). 

What remains to be calculated is the correc- 
tion coming from real emissions, confined to the 
band (approximately) between the shaded and 
dotted lines. As long as the observable is such 
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that one can ensure that this band has a width 
lnl/e that is not parametrically large (i.e. inde- 
pendent of L), then the band will give at most 
single- logarithmic contributions a"L n , i.e. it will 
contribute only to the g2(a a L) function, not to 
g\. Because the density of primary emissions per 
unit rapidity in the band, ~ ot s lnl/e will be low, 
emissions will be widely separated in rapidity, al- 
lowing one to make use of the property of an- 
gular ordering whereby the full matrix element 
for multiple emission reduces to that for indepen- 
dent emission 0, simplifying the calculation of 
gi- Note that this only works for continuously 
global observables |4I5| . 

To aid the mathematical treatment of the 
above procedures it is convenient to parametrise 
the observable's dependence on a single emission 
close to hard leg £ (in e + e _ — > 2 jets, £ = 1,2) as 



V(k)=d m (4,)(^-) e 



(5) 



with k t and r\ now measured with respect to leg 
£, and a, bt and de numerical constants and ge(<fr) 
a function of the azimuthal angle 0, all of them 
characteristic of the observable and relevant to 
determining the position of the dashed line in figH 
for a given v (where <f> is not represented). One 
also defines dg = de exp [J^ 71 |^ In <#($)]. 

It will be useful to express momenta in terms 
of the effect they have on the observable, through 
functions that satisfy 



(6) 



This condition alone is not enough to fully specify 
«i(C) — additionally one needs to know the leg £i 
to which it is collinear, its azimuthal angle and 
how its rapidity depends on £ — for concreteness 
here we will take this to be (C) = & (In. 1/C) / 
be) so that taking & in the range < & < 1, «i(C) 
will span from the large-angle soft region to the 
hard-collinear region (regardless of £). 

Part of the logic of this scaling is that the emis- 
sion pattern ~ a s ^dr/ becomes proportional to 
a s L ■ dln£<i£, with the single logarithmic factor 
a s L factoring out explicitly. It also helps com- 
pactify conditions such as eq. which becomes 



V(ki((i), . . . ,«n(Cn)) ~ max(Ci, ■ • - C«) • 



(7) 



Actually "~" is difficult to encode directly in a 
computer, as will be needed for automation. Here 
the meaning of "~" is that the ratio between the 
two sides should not be parametrically large, even 
if one scales all the Q by a common large fac- 
tor (i.e. no new large ratios should appear in the 
problem, because they could lead to large loga- 
rithms). This suggests the following more precise 
condition, 



lim 



V(«i(uCl)) ■ ■ M«n(«Cn)) 



— /{k«}(Ci> ■ ■ ■ , Cre) ) 

(8) 



where f{ Ki } should be some well-defined non-zero 
function and where the subscript indicates that it 
may depend also on the additional characteristics 
of the Ki (£i, £i, ef>i). 

The also useful in formulating the pre- 

cise conditions needed to go from figCJi to figGJx 
For example we needed the observable to be such 
that we could neglect the grey emissions, leaving 
only those in the band between the dashed and 
dotted lines — furthermore we needed to be able 
to require that the width of the band, In 1/e, was 
not parametrically large, i.e. would not be associ- 
ated with extra logarithms. This is equivalent to 
saying that in eq. (JSJ), independently of v on the 
left-hand-side, there should exist an e such that 
any emissions with Q < e can be neglected, or 

f 1 ™ Q /{«<}(&> •••iCn) = /{ K< }(Cl: ■ ■ -jCn-l) ■ (9) 

This looks very much like infrared safety, but 
is stronger because /{ Ki } already involves an in- 
frared limit of its own. Accordingly we call it 
recursive infrared safety. Making use simulta- 
neously of normal infrared safety it can also be 
rewritten in terms of a commutator of limits, 



lim , lim 

C„^0 v^O 



V(Kl(vCl), K n (v( n )) 



= 0. 



(10) 



Similarly, to ensure that one can replace the clus- 
ters of black emission in fig^i with single emis- 
sions in fig^D, one needs a condition called recur- 
sive collinear safety. 
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Given these recursive IRC (rlRC) conditions, 
exponentiation is guaranteed, and to NLL accu- 
racy one has 

lnSy(w) = Lg 1 (a s L,a,{b i })+ 

+ g 2 {a s L, a, {be}, {dt}) + T V {R') , (11) 

where the Qi are known analytical func- 
tions whose observable dependence enters only 
through the a, be and de coefficients; R' = 
d L {LQi{a s L,a,{b t })) and the NLL T V {R') ac- 
counts for the observable's non-trivial dependence 
on multiple emissions. In most cases it is given 
by 

R' 00 1 

F v (R') = ]im— x 

m=0 

(m+l n .1 .» r 2ir ji \ 

nix/fi £)w 

x exp(-i?'ln/ {K!} (Ci,...,Cm+i)) ■ (12) 

where R' e . is the contribution to R' from leg I. 

Resummation is automated through an expert 
system (CAESAR jSj) which uses high precision 
arithmetic to (a) probe the observable in the soft 
and collinear region to obtain the a, be and de 
coefficients and the ge(4>), (b) test the rlRC safety, 
(c) evaluate J-'v(R') via Monte Carlo integration. 

One of the most non-trivial applications of 
CAESAR is to dijet event shapes at proton- 
(anti)proton colliders. Such event shapes are of 
interest for various reasons, among them sensitiv- 
ity to a new class of soft colour-evolution anoma- 
lous dimensions [7], and the prospect of extract- 
ing information about non-pcrturbative correc- 
tion from gluon hadronisation and from the un- 
derlying event. Unfortunately the event shapes 
measured so far [5] are not global and therefore 
not within the scope of CAESAR. Globalness is 
a significant issue at hadron-colliders because de- 
tectors are not able to measure the whole event, 
but only up to some limited maximum rapidity, 
?7 max ~ 3.5 at the Tevatron (~ 5 for LHC). 

There are two main ways of working around 
this problem. Firstly one can define an event 
shape that measures everywhere (directly global) 



[5], and then note that if ?7 max is sufficiently large, 
then the error introduced by actually measur- 
ing only particles with \n\ < ?7 max is of order 
SV = d i e - ( a+f ' i )' 7max , where bj, d~i are the be, de 
coefficients for the incoming hadron legs. As long 
as most of the event shape distribution is con- 
centrated at values of the observable much larger 
than this, then the non-measurement of particles 
with 1 77 1 > ?7 max has a negligible impact on the 
distribution. 




ln(T mg ) 

Figure 2. Resummed distribution of T ms for 
Tevatron events where the hardest jet has E± > 
50GeV. The point where the 7? max = 3.5 cutoff 
becomes relevant is indicated by a hashed line. 

Among these directly-global observables, there 
are two subclasses. Firstly there are those nat- 
urally defined in the whole event, a transverse 
thrust t_l )9 , its minor T m ^ g and the resolution pa- 
rameter at which the k t jet-finder goes from clas- 
sifying an event as 2 + 2-jet like to 3 + 2-jet like, 
2/23- The resummed distribution for T mj9 is shown 
in fig. 

Alternatively one can define a main event 
shape in a central region C {e.g. \n\ < 1) and 
add to it an exponentially-suppressed contribu- 
tion for particles in the forward region £q = 
J2 l ^c k ti/Q±,ce^ m " ,,c K where Q ±fi is the sum 
of transverse momenta of particles in C and r\c 
is the (transverse-momentum weighted) mean ra- 
pidity of particles in C. Such 'combined' observ- 
ables can be constructed starting from the 'nat- 
urally global' event shapes given above, and also 
from event shapes that only make sense when de- 
fined in the central region such as jet broadenings 



Automated Resummation and Hadron Collider Event Shapes 



5 



Table 1 

Properties of the various event shapes proposed in |10j . Entries marked with an asterisk are subject to 
uncertainties. 



Event-shape 


Impact of ?7 max 


Resummation 
breakdown 


Underlying 
Event 


Jet hadronisation 




tolerable* 


none 




-i/Q 


T m ,g 


tolerable 


none 




~ 1/(V57Q) 


2/23 


tolerable 


none 


~ VV23/Q* 


~ v^/Q* 


T±,e, px,e 


negligible 


none 


-i/Q 




Bx,s 


negligible 


none 


-l/Q 


- 1/(VS7Q) 


T m ,£ 


negligible 


serious 


-l/Q 


- i/(ys7Q) 


V23X 


negligible 


none 


-1/Q 




t_l,ii, px,n 


none 


serious 




-i/Q 


T m ,TZ, Bx,TZ 


none 


tolerable 


~1/Q 


~ 1/(VS7Q) 


y '23, 


none 


intermediate* 


~ V2/23/Q* 





and invariant masses, giving e.g. Btx- The nat- 
urally global observables always have bi = and 
are more sensitive to the ?7 max cut (cf. fig.0) than 
those combined with the forward exponential sup- 
pression, which have b t = a. 

One can also define 'indirectly-global' observ- 
ables. Here too there is a main event shape in the 
central region and to it one adds a recoil term 



/Q±,C measured in C, which 
through momentum-conservation is sensitive to 
emissions outside C. This eliminates the problem 
of sensitivity to ?7 max , however instead one runs 
into the issue that the CAESAR resummation is 
reliable only as long as the main mechanism by 
which an observable can take a small value is by 
suppression of radiation. If instead it can have a 
small value through cancellations between differ- 
ent emissions (as with a vector recoil) then the re- 
summation breaks down via a divergence oiJ-(R') 
at some small but finite value of v. 

Table ^ summarises the impact of 77 max and 
resummation breakdown (if any) for a range of 
observables defined in ^U]. It is gives expecta- 
tions for the degree of sensitivity to the underly- 
ing event and jet hadronisation, illustrating the 
complementarity between observables. 

The work presented here was carried out in col- 
laboration with A. Banfi and G. Zandcrighi. I 
wish to thank the organisers of the RADCOR 
2005 Symposium for the welcoming and stimu- 



lating atmosphere of the conference as well as for 
financial support. 
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